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We study the effect of disorder on massless, spinful Dirac fermions in two spatial dimensions 
with attractive interactions, and show that the combination of disorder and attractive interactions 
is deadly to the Dirac semimetal phase. First, we derive the zero temperature phase diagram of 
a clean Dirac fermion system with tunable doping level (fi) and attraction strength (<?). We show 
that it contains two phases: a superconductor and a Dirac semimetal. Then, we add disorder, and 
show that arbitrarily weak disorder destroys the Dirac semimetal, turning it into a superconductor 
instead. Thus for Dirac fermions near charge neutrality, disorder actually assists superconductivity. 
We discuss the strength of the superconductivity for both long range and short range disorder. 
For long range disorder, the superconductivity is exponentially weak in the disorder strength. For 
short range disorder, a uniform mean field analysis predicts that superconductivity should be doubly 
exponentially weak in the disorder strength. However, a more careful treatment of mesoscopic fluc- 
tuations suggests that locally superconducting puddles should form at a much higher temperature, 
and should establish global phase coherence at a temperature that is only exponentially small in 
weak disorder. Thus, mesoscopic fluctuations exponentially enhance the superconducting critical 
temperature. We also discuss the effect of disorder on the quantum critical point of the clean sys- 
tem, building in the effect of disorder through a replica field theory. We show that disorder is a 
relevant perturbation to the supersymmetric quantum critical point. We expect that in the pres- 
ence of attractive interactions, the flow away from the critical point ends up in the superconducting 
phase, although flrm conclusions cannot be drawn since the renormalization group analysis flows to 
strong coupling. We argue that although we expect the quantum critical point to get buried under 
a superconducting phase, signatures of the critical point may be visible in the flnite temperature 
quantum critical regime. Our results have implications for experiments on proximity induced su- 
perconductivity in Dirac fermion systems, where they imply an enormous disorder-enhancement of 
the superconducting susceptibility. As a result, the proximity induced superconductivity in dirty 
systems is expected to be much stronger than that in clean systems at the Dirac point. 



I. INTRODUCTION 

The study of many body effects in disordered elec- 
tronic systems has been a focus of condensed matter re- 
search for decades.'i^ An important subset of problems 
involves disordered electronic systems with attractive in- 
teractions, and hence the interplay of superconductivity 
and disorder.-^ The Anderson theorenP states that s-wave 
spin singlet superconductivity is immune to weak time 
reversal invariant disorder. However, disorder is widely 
believed to suppress all other forms of superconductivity. 
The Abrikosov-Gorkov theorjE' provides a framework for 
analyzing the suppression of superconductivity by disor- 
der. More recently, it has been pointed out that meso- 
scopic fluctuation effect^^^^^ can render superconductiv- 
ity more robust to weak disorder than the Abrikosov- 
Gorkov theory would predict, but disorder nevertheless 
suppresses superconductivity. 

The recent discovery of two-dimensional (2D) Dirac 
fermion systems, such as graphen^ or the surface of a 
3D topological insulator,'^ has opened a new avenue of re- 
search into the interplay of disorder and superconductiv- 
ity. The study of disordered Dirac fermions began almost 
30 years agcPE] and it is believed that a single species 
of Dirac fermions is protected against localization,* al- 
though a random scalar p otent ial disorder generates a 
non-zero density of states.f^iEI!, Meanwhile, it is also 



known that Dirac fermion systems at charge neutrality 
do not develop superconductivity for arbitrarily weak at- 
tractive interactions, and that there is a quantum criti- 
cal point at a critical attraction strength which separates 
the Dirac semimetal from the superconductor^^t^ In a 
recent stimulating development, it has been pointed out 
that this quantum critical point has an interesting ef- 
fective field theory descrip tion, which displays emergent 
dynamical supersymmetry.^iSlin! However, while disorder 
and attractive interactions have been studied in isolation 
for topological insulator surface states, the combination 
of disorder and attractive interactions has not been stud- 
ied. 

In this paper, we study Dirac fermions with attrac- 
tive interactions, in the presence of disorder. We show 
that the combination of scalar potential disorder and at- 
tractive interactions is particularly deadly to the Dirac 
semimetal, driving a transition into a superconducting 
phase for arbitrarily weak attraction strengths. Remark- 
ably, for Dirac fermion systems disorder actually en- 
hances superconductivity, allowing superconductivity to 
develop where the clean system would have been semi 
metallic. We also show that while the disorder enhance- 
ment of superconductivity can be estimated using a mean 
field theory a la Abrikosov-Gorkov, this treatment dra- 
matically underestimates the strength of the supercon- 
ductivity. A proper treatment of mesoscopic fluctuations 
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is necessary to determine the strength of superconduc- 
tivity in the disordered Dirac fermion system. Our re- 
sults may also have implications for Dirac fermion sys- 
tems with repulsive interactions and disordei^^. Dirac 
fermion systems with repulsive interactions are described 
by a Gross Neveu model, which also has a quantum criti- 
cal pointPSl which may be destabilized by disorder. How- 
ever, we do not pursue this line of research in this pa- 
per, leaving it as a topic for further work. We note too 
that whereas this paper is focused on Dirac semimetals, 
the interplay of disorder and superconductivity in normal 
metals has been studied in Ref. [23J 

While Dirac fermions arise both on the surface of a 
topological insulator and in graphene, in graphene the 
Dirac fermions are fourfold degenerate, coming in two 
spin and two valley flavors. On the surface of a topolog- 
ical insulator, however, there is just a single species of 
Dirac fermion. In this paper, we focus our analysis on 
the case of a topological insulator, with a single species 
of Dirac fermion. This captures the essential physics 
of Dirac fermions with attractive interactions and dis- 
order, but is easier to treat analytically because of the 
smaller number of degrees of freedom. The basic con- 
clusions should also apply to graphene, insofar as disor- 
der will generate a non-vanishing density of states and 
enable superconductivity to develop for arbitrarily weak 
interactions. However, the existence of a valley degree 
of freedom in graphene, and the fact that disorder can 
cause intervalley scattering, may lead to additional fea- 
tures not present in the problem studied here. We leave 
the generalization of this analysis to graphene as a topic 
for future work. 

In this work, we consider a model of Dirac 
fermions with purely attractive, phonon mediated inter- 
actions, neglecting the Coulomb repulsion. In princi- 
ple, the (unscreened) Coulomb interaction can prevent 
superconductivitji2l^ and the study of superconductiv- 
ity in the presence of Coulomb interactions is a sub- 
ject we leave for future work. However, we note that 
in experiments on graphene or topological insulators, the 
Coulomb interaction can always be screened by metallic 
gates, so the neglect of the Coulomb repulsion can be 
justified. 

We also emphasize that our work has important im- 
plications for ongoing experiments attempting to induce 
superconductivity in graphene and topological insulators 
by means of the proximity effect. In the context of the 
proximity effect, the disorder enhancement of Tc which 
we identify can be read as a disorder enhancement of the 
superconducting susceptibility. Thus, even in materials 
where there may not be an intrinsic attraction or intrin- 
sic superconductivity, the proximity induced supercon- 
ductivity in the disordered system will be dramatically 
stronger than than in the clean system, for reasons we 
explain in this paper. 

This paper is structured as follows. In Sec. II, 
we discuss the zero temperature phase diagram of the 
clean Dirac fermion system in the presence of attrac- 



tive delta function interactions. While the phase struc- 
ture at the Dirac point (chemical potential /z = 0) as 
a function of attra ction strength g has previously been 
understood, 'i^E^'2M28l|;]2g phase diagram in the 5 — /i 
plane has not been presented as far as we are aware. In- 
deed, the — >■ cxD limit remains controversial, with the 
existing literatur^^MSZJ apparent disagreement. We re- 
solve this disagreement by means of a careful analysis 
that takes into account the finite ultraviolet cutoff for 
the interaction, which justifies use of a projected Hamil- 
tonian. The projection operation introduces a gauge re- 
dundancy, which must be dealt with carefully. We also 
show how the system interpolates between the = 
and large /j, limits, and present the full zero temperature 
phase diagram. The zero temperature phase diagram 
contains two phases: a Dirac semimetal along the fi = 0, 
g < gc line segment, and a superconductor everywhere 
else. There is a single superconducting phase, which is 
fully gapped and preserves time reversal symmetry every- 
where. For small doping, the superconductivity is mostly 
spin singlet and s-wave. 

In Sec. HI, we discuss the influence of disorder. First, 
we discuss disorder in the attraction strength g. We show 
that this form of disorder is a relevant perturbation at the 
fi — 0, g — gc critical point, and must necessarily change 
the universality class. However, we expect disorder in 
the attraction strength to be weak at the bare level, and 
thus the effect of attraction strength disorder may not 
manifest itself until very low energy scales. However, 
in any realistic experimental sample, there is likely to 
be significant disorder in the chemical potential /z. We 
show that smooth chemical potential disorder destroys 
the semimetal phase, producing a network of electron and 
hole doped puddles. In the presence of attractive inter- 
actions, this system develops percolating superconduc- 
tivity for arbitrarily weak attractive interactions, with 
a critical temperature that is exponentially small in the 
typical doping. We also show that short range disorder 
(which we model in terms of delta function impurities) 
also destroys the semimetal phase, introducing a low en- 
ergy density of states that is exponentially small in the 
disorder strength. A uniform mean field analysis of the 
pairing problem suggests that the short range disordered 
Dirac fermion system should display superconductivity 
for arbitrarily weak interactions, but with a critical tem- 
perature that is doubly exponentially weak in the disor- 
der. 

In Sec. IV, we analyze the influence of mesoscopic fluc- 
tuations on superconductivity for short range disorder, 
in the weak coupling limit. Our analysis suggests that 
mesoscopic fluctuations dramatically enhance supercon- 
ductivity. In particular, the analysis suggests that small 
puddles of local superconductivity appear at tempera- 
tures that are linearly small in weak disorder, and the 
Josephson coupling between these puddles establishes 
global phase coherence at a temperature that is expo- 
nentially small in weak disorder. This represents a dra- 
matic enhancement over the uniform mean fleld theory. 
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which predicts a critical temperature that is doubly ex- 
ponentially small in weak disorder. The application of 
a transverse magnetic field frustrates the global phase 
coherence, driving a phase transition into a gauge glass 
phase. 

In Sec. V, we analyze the interplay of disorder and at- 
traction within a rcnormalization group (RG) framework, 
in the vicinity of the quantum critical point ^ — Q,g — Qc- 
The construction of an RG for the attractive interaction 
requires us to work in an e-expansion about the upper 
critical spacetime dimension D = A. Meanwile, disorder 
is taken into account through a replica field theory ap- 
proach. The discussion in this section complements the 
discussion in Sec. Ill and IV. Whereas Sec. Ill discussed 
the interplay of disorder and interactions near the free 
fermion point, the renormalization group analysis stud- 
ies the interplay of disorder and interactions near the 
quantum critical point. We find that whereas chemical 
potential disorder is a marginally relevant perturbation 
at the Gaussian point, it is a (power law) irrelevant per- 
turbation at the quantum critical point. Thus, one might 
naively expect the quantum critical point to be stable in 
the presence of weak chemical potential disorder. How- 
ever, a careful analysis reveals that in fact chemical po- 
tential disorder is a dangerous irrelevant perturbation, in 
that it generates disorder in the BCS coupling, which is a 
relevant perturbation to the critical point. Thus, the su- 
persymmetric critical point of Ref . [TS 20 is in fact unsta- 
ble in the presence of disorder. We expect that the flow 
away from the critical point leads to the superconducting 
phase, although firm conclusions cannot be drawn since 
the RG flows to strong coupling. 

Finally, in Sec. VI we discuss the prospects of observing 
signatures of the supersymmetric quantum critical point 
identified in Ref. 1191201 given the inevitability of disor- 
der. We point out that even though the quantum critical 
point will be buried under the superconducting phase in 
any realistic experimental sample, in a sufficiently clean 
sample some signatures of the critical point may survive 
in the finite temperature quantum critical regime. We 
conclude by discussing future directions for the study of 
disorder and attractive interactions in Dirac fermion sys- 
tems. 

We note that previous studies of Dirac quasiparticles 
in nodal superconductors have found a (secondary) su- 
perconducting transition in the disordered system, when 
no such transition occurs in the clean one^SHSSl. There 
is some mathematical resemblance between these results 
and the discussion in Section III. However, the colossal 
enhancement of superconductivity by rare region effects 
discussed in Section IV, as well as the strong coupling 
physics discussed in Section V,VI, have no analog in these 
works. The discussion in section HI also has important 
differences, in that it belongs to a different Cartan sym- 
metry clasJ^, with very different localization physics, 
and also in that we are discussing a true superconducting 
instability in a semimetal, whereas the work^SlHlll g^j-g (jjg. 
cussing the appearance of a secondary s-wave component 



of the order parameter in a d-wave superconductor. 

II. PHASE DIAGRAM OF CLEAN DIRAC 
FERMIONS WITH ATTRACTIVE 
INTERACTIONS 

We begin by considering a single species of Dirac 
fermions in the absence of disorder, with short range at- 
tractive interactions. The creation operator for Dirac 
fermions with momentum /ik is V'k- The spinor struc- 
ture of the Dirac fermions is implemented by defining 
V'k ~ ("^t k' ''1 k)' where ^ creates a fermion with spin 
a and wavevector k. It is convenient to introduce the 
Pauli matrices Gi which act in the spin space, and to 
also define ctq to be the two dimensional identity matrix 
acting in spin space. The second quantized Hamiltonian 
may then be written as H = Hq + Hi , where 

Ho = ^ipi{-fJ.cro + vkx(Ti+vkya2)ipu, (1) 
k 

= ^(q'l^'P)V'Icro'0k+qV'J,O'oV'p-q, (2) 

k,p,q 

where v is the Fermi velocity, fi is the chemical potential 
(which controls the doping level), and we have assumed 
a purely 'density-density' interaction, which has no in- 
trinsic spin structure. The interaction is assumed to be 
attractive, V < 0. Although we have taken ct to be a 
spin index for simplicity, we are aware that for generic 
topological insulator surfaces it may be a composite spin- 
sublattice indejP^. The distinction makes no difference 
to our analysis, which involves pairing between time re- 
versed states (not necessarily opposite spin states). The 
Pauli matrices should thus be understood as acting in 
the spin/pseudospin space relevant for the surface states, 
such that a and —a are time reversed states. 

We now project onto the Cooper channel, by restricting 
the interaction Hamiltonian Q to p = — k. The projec- 
tion is appropriate for studying the superconductivity of 
Dirac fermions. The projection leads to a Hamiltonian 
of the form Ha + H[ , where Hq is given by ([T]) and 

H'l=Y. ^(q' k)^I'To^qVlk^O^-q- (3) 

k,q 

We further restrict ourselves to a short range, 5 function 
interaction. In this case, V is independent of momenta 
and we obtain the BCS Hamiltonian for Dirac fermions 

~ sJI^k^o^q^-k^oV'-q, (4) 

k,q 

where g = is the superconducting coupling. The rest 
of this section will be devoted to solving Q. 

We proceed as follows. First, we provide a general dis- 
cussion aimed at classifying potential solutions of m into 
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distinct phases. Then we solve Q along the line n = 0. 
The solution along this line is well known, but we present 
it for completeness sake. Next, we discuss superconduc- 
tivity in the limit \fi\ — ^ oo. Our discussion resolves a 
disagreement between Ref. and Rcf. 26.271 as to the 
nature of superconductivity in this limit, and highlights 
subtleties connected with UV cutoffs and gauge ambigu- 
ities which must be properly taken into account to un- 
derstand this limit. Finally, we solve Q for arbitrary 
doping, and construct the full phase diagram in /i — .g 
space. 

We show that the entire fi — g plane has a ground state 
that belongs to a single superconducting phase, with the 
exception of a single line segment along the /i = line, 
which is semi-metallic. This follows because the Hamil- 
tonian Q has a superconducting ground state for any 
system with a non- vanishing density of states, for arbi- 
trarily weak interactions, and the low energy density of 
states vanishes only on the fi — line. More non-trivial 
is the fact that the /i > and fi < superconducting 
regions belong in the same phase - a phase which is fully 
gapped and time reversal invariant, with an order pa- 
rameter that is a real linear combination of spin singlet 
and spin triplet parts. A spin triplet component emerges 
because the spin singlet order induces spin triplet order 
at any non-zero doping /i 7^ 0. 

A. Symmetries and superconductivity 

The most general order parameter contains spin singlet 
and triplet pieces, and can be written as 

k 

Fermi statistics demand that A^^k = As^_k and dk = 
— d_k- The Hamiltonian is invariant under time reversal 
(TRS), which is implemented by complex conjugation, 
taking k — >■ — k, and acting with the operator ia2. Thus, 
we conclude that TRS is preserved if and only if dk/A^ 
is a vector with purely real components (i.e. if the super- 
conductivity is a real linear combination of singlet and 
triplet) . There is also a particle hole symmetry along the 
/i = line and we will discuss this symmetry when ap- 
propriate. However, this symmetry is broken by scalar 
potential disorder (our primary focus in this paper), and 
so we do not attach much weight to this symmetry. How- 
ever, there is a second mean ingful distinction between 
superconducting phases'^^"'^^' namely whether there ex- 
ist gapless Bogoliubov-de Gennes (BdG) quasiparticles. 
As we shall see, the superconducting solutions to Q are 
always fully gapped and invariant under time reversal, 
so the phase diagram contains a single superconducting 
phase. 

We assume that spin singlet superconductivity is the 
primary superconducting instability. We now discuss 
whether spin singlet order can induce spin triplet order. 
At non-zero fj,, the only symmetries are TRS (discussed 



above), and also a continuous rotation symmetry imple- 
mented by the generator = Lz + ^a^. Spin singlet 
order is even under TRS, and has Jz — 0- Time reversal 
symmetry prohibits the spin singlet order from coupling 
to any spin triplet order parameter with complex dk/Aj, 
whereas rotation symmetry prohibits it from coupling to 
any spin triplet order with Jz 7^ 0. However, a spin 
triplet order parameter with dk/A^ = k is even under 
TRS, has Jz = 0, and can thus couple directly to the 
spin singlet order parameter. Thus, in general we expect 
a solution of the form 

A = ^(V'k(A, -KF(^)k-a)ia2V^_k), (5) 

k 

where k = [k^, ky) and the spin triplet piece of the or- 
der parameter is induced by the spin singlet piece. The 
proportionality constant F{^) remains to be determined. 

The Dirac point /i = is special, in that a large number 
of extra symmetries appear. In particular, the theory 
becomes Lorentz invariant, and also there exists at /i = 
a particle hole symmetry,!^ under which spin singlet and 
spin triplet pairing are even and odd respectively. Thus, 
the spin singlet and spin triplet orders cannot couple at 
/I = 0, and it follows that F(/x = 0) = 0. 

The remainder of this section will be devoted to ex- 
plicitly obtaining and solving BdG equations leading to 
Eq.([5]). We will show that F{fj, — > ±00) — sign(/x) and 
F(/i — >■ 0) ~ g'fi, where g' is the (/i dependent) attrac- 
tion in the triplet channel. Readers uninterested in the 
clean system details may skip directly to Sec. IH. 

B. Superconductivity along the ^ = line 

We wish to solve H^cs at /i = 0. Since the Pauli 
principle ensures that two fermions with the same spin 
cannot interact through a delta-function potential, we 
can re-write the Hamiltonian as 

H = il^l f ^fcx^i + vkya2 ] tpk 

k ^ ^ 

- 5X!'^k(-*'^2)V'lk'/'q«CT2V'-q- (6) 
k,q 

Superconductivity is necessarily spin singlet, and is char- 
acterized by the order parameter 

A. -5](V^k(^a2)V;-k). (7) 

k 

It is now convenient to introduce the Euclidean time 
path integral representation of the partition function, 
Z = J i:)[V'^V']exp( - /jf dr / d^xL[iP\ij]), where /3 is 
the inverse temperature and the Lagrangian takes the 
form 
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where we have replaced the sum over momenta by 
an integral. Since we are working in the path in- 
tegral representation, tp^ and ijj now represent Grass- 
man valued fields rather than second quantized op- 
erators. We define the four component Nambu 
spinor 4* = (V'k, V'lk)- After decoupling the four- 
fermion interaction in ^ by means of a Hubbard- 
Stratonovich transformation, the partition function 
can be rewritten as Z = J D[^\^ , A* , A\exp [ - 

Jq dr J d'^xL[i>'^,'ii,A*,A]), where A is a complex val- 
ued (bosonic) field, and 



L = T 



(27r) 



I A. 



G. 



/ 




vk^ 









iuin - M 


-A, 


I 





-A* 


iUJn + 


a: 





vk^ 



2<? ' 

A, 




where we have introduced the Matsubara frequencies 
ujn — (2n + l)7rr, and have assumed that the s-wave 
order parameter A^ is isotropic in momentum space. 

We can now integrate out fermions exactly to obtain a 
Lagrangian that only involves the order parameter fields. 
This Lagrangian takes the form 



L = 



d^k 

d^k 



Tr In _k 
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ln(c.2+«2^.2 + |A,n2 + 



I A, 



2<? 



Variation with respect to |Ap then yields the gap equa- 
tion, 



2T 
2^ 



E 



kdk 



1 



The momentum integrals carry an implicit cutoff at the 
scale vk = wu, where w^i is the Debye frequency. At zero 
temperature 
be solved as 



J 1^ , and the gap equation can 



which has solutions only for 

g> 9c = 2t:v^/ujd- 

Thus we recover the well known result that supercon- 
ductivity for Dirac fermions along the /i = line is a 
threshold phenomenon, with the order parameter devel- 
oping a non-zero expectation value only if interactions 
are strong enough. The solutions (assuming |As| <^ujd) 
take the form 



2'KV^ 



2ttv^ , 



where the approximate equality holds close to the thresh- 
old, 5g = g — gc gc, but not so close as to be governed 
by the critical point. Thus, we conclude that along the 
fi — line, the system has spin singlet superconducting 
order for g > gc and is a semimetal for g < gc- As has 
been pointed out in Rcf. I19|201 the critical point g = gc 
has some unusual features, and is described by an effec- 
tive field theory that exhibits emergent dynamical super- 
symmetry. 



C. Superconductivity at large doping, ujd 

At large doping, the spin basis is not the most conve- 
nient basis to work with, since the low energy states are 
linear superpositions of spin up and spin down states. In- 
stead, we transform to a basis of + and — helicity states 
(upper and lower Dirac cones) by performing a unitary 
transformation, according to 



U 



c+,k 
c-,k 



V2 



_g-i0k/2 



(8) 



where k^ + iky — 



lj|g«0k^ The requirement that U should 
be single valued demands that G should be a half inte- 
ger. We wish to stress that the choice of unitary matrix 
in ([s]) is not unique, since we can freely choose G to be 
any half integer. This 'gauge ambiguity' makes no differ- 
ence if we work with the full Hamiltonian. However, we 
will shortly be projecting onto a single helicity basis, and 
the projected Hamiltonian will look different with differ- 
ent gauge choices. The different projected Hamiltonians 
should be understood as being gauge equivalent. 

We now use the unitary transformation (|8| to express 
the Hamiltonian in the helicity basis. Defining tq to 
be the identity matrix in helicity space, and Ti to be 
the Pauli matrices in helicity space, we can rewrite the 
Hamiltonian in the helicity basis as 



k k,q 



,2iG(0q-0k) 



X $ 
X $ 



t 
k 

t 

-k 



2 cos 



2 cos 



To 4- 2^ sm — - ti 



^k — 0q \ , o ■ ■ / V^k ~ V'q 

' To + 2i sm ' 



q' 
(9) 



where we have not yet specified the choice of gauge G. 

It is intuitively obvious that at large doping, only 
states close to the Fermi surface need to be considered, 
and thus one can project onto the states with helicity 
sign(/i). To justify this projection, we note that the in- 
teraction implicitly has an ultraviolet cutoff on the scale 
ujd. Thus, it cannot couple states near the Fermi surface 
to states further than ujo away from the Fermi surface. 
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If the doping > u;^), then one can project onto states 
with a single hehcity. The apparent neglect of this UV 
cutoff in Ref. [55] explains the discrepancy between the 
large doping results in Ref. [55] and Ref. I26I27I 

We consider electron doping ^ > 0. The case of hole 
doping /i < follows by analogy. After projection onto 
states with positive helicity, the Hamiltonian becomes 



H 



k 



.( — ^ + vk)c+,] 



k,q 



^ '^+,k'^+.q'^+,-k^+.-q- 

Now cos2a;/2 = i(l + cos(a;)) = i(2 + + e"*^). Thus, 
note that the projected interaction has harmonics with 
angular momentum 2G, 2G + 1, 2G — 1, where G can be 
any half integer. It is most convenient to make the gauge 
choice G = 1/2. Then the attractive potential has har- 
monics with angular momenta I = 0,1,2. Note that an 
effective 'p-wave' (1=1) harmonic has been generated by 
the projection, even though we started with a purely s- 
wave interaction.l^Sl 

Now since after projection we are dealing with a stan- 
dard one band BCS problem for spinless fermions, Fermi 
statistics demand that the superconductivity has to be 
odd parity i.e. A_k = ^ Ak- Thus, the / = 0, 2 har- 
monics do not introduce superconductivity, and may be 
projected out. We need retain only the p = 1 harmonic, 
which gives us a HamiltoniarPSl 



H 



-X 



k 

k,q 



kC+,qC_|_ _k'^+,-q- 



Proceeding to the path integral representation, and de- 
coupling the four ferniion interaction using an order pa- 



rameter field A_|_ — I J2q 
clidean time Lagrangian 



;_n)e*'*'i, we obtain the Eu- 



k 

-f ^(A+CqC_qe'^o+c.c.) + i|A+p. 



and summing over fermonic Matsubara frequencies gives 



(2^)2 2^/{^ik^^ 



^^^= tanh I — 



where the integral goes over fi — ujd < vk < ii + ujjj, and 
we are working in the limit fi 3> ujd- In this limit, the 
zero temperature gap equation can be solved as 



In 



27rw2 
.9M 



In the weak coupling limit 2ttv'^ jg^, ^ 1, this gives us 



|A+| = 2w_DCxp - 



(10) 



Thus one concludes that there is an order parameter 
A+(k) = I Sn(cqC-q)e"^'' with expectation value given 



by ( 10 ). For a different gauge choice, the expression ( 10 ) 



still holds, but the order parameter instead has the form 
A+(k) = i Eq(cqC-q)e*^<=^'^''. The calculation for < 
is analogous, with the only difference being that the pro- 
jection is now on the band with negative helicity. Thus, 
we conclude that in the helicity basis, the order parame- 
ter takes the form 



A = ^ e^'"^*^ 



($J(ro -f sign(/i)T3)$_ 



The BdG spectrum is fully gapped. Meanwhile, trans- 
forming back to the spin basis, the order parameter be- 
comes 



sign(/i)(cos ^iikO"! - sino-2))«cr2*-k) 



(11) 

Note that this is a real linear combination of spin singlet 
and spin triplet terms, and thus preserves time rever- 
sal symmetry. Thus, the superconducting state at large 
doping is in the same phase (fully gapped, time rever- 
sal invariant) as the superconducting phase at ^ = 0. It 
seems likely that the system should smoothly interpolate 
between the /x = and large /i limits. We confirm this 
in the next section. 
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After integrating out the fermions, we obtain an action 
purely in terms of order parameter fields, which takes the 
form 



L = -Tr In 



Lol + {vk-^if + \A+\' 



1 



where Tr denotes summation over Matsubara frequencies 
and integration over momenta. Variation with respect to 
|A_|_p yields the gap equation, 

Tr 



Lol + {vk-^JiY + \^+\^ g' 



D. Superconductivity at small, nonzero ^ 

We have solved the Hamiltonian Q in the fi — and 
large fi limits. We now turn to the small fi regime. We 
have seen that a finite UV cutoff ujjj on the interac- 
tion generates a spin triplet component to the interaction 
away from — 0. We therefore introduce a weak spin 
triplet component into the action by hand, and write the 
Lagrangian (after Hubbard-Stratanovich decomposition) 
as 
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iujn — fJ. vk^ 

vk^ iuJn — fJ. 

-dlik)-zdUk) -A*+dUk) 

A*,+dl{k) dl{k)-id*2(k) 



-rfi(k)+zd2(k) A,+d3(k) 

-As+d3{k) di(k)+zd2(k) 

iiOn + fJ- vk^ 

vk+ iuJn + /i 



2ff 



25' 



r 



where g' is a function of Consistency with the 

/i = and large /i Hmits requires that g'(0) = and 
g'{\^J.\ > ujd) — g- We assume that g' interpolates 
smoothly in between, so that < g' < g in the small 



/i regime. Meanwhile, the Matsubara frequencies are im- 
plicitly summed over, and di,2,3 are odd functions of k. 

After integrating out fermions, we obtain 



L = -Tr In ^ ^2^2 ^ ^^^^2 ^ fc2)|^^|2 ^ 2(c^2 _ A:2)|d3(k)|2 + 2(c^2 _ ^2 ^ fc^)|^2(]^)|2 ^ 2(c^2 ^ ^2 _ fc2)|^^(j^)|2 

3 

-4fc^fcj,(di(k)d;(k) + d*(k)d2(k)) + Awk^{d^{k)dl{k) + c.c.) + 4.wky(dz{\<i)d\{k) + c.c.) + lA,]"* + ^ M^(k)|^ 



A^[K(k)f + K(k))2 + (d5(k))2]+c.c. 



di(k)^4(k)' + di(k)'d;(k)2 + d3(k)2d;(k)2 + c.c 



~4fcj,A*(A:d2(k) + A,d;(k)) + 4A:,Ai(A:di(k) + A,rf*i(k)) + 2ixWAs\^ + Y,dj + lo"^ - k") + // 



2i7 



25' 



(12) 



where we have set v = 1 for simplicity. Now, g > g' in 
this regime, so that Ag is the primary instability. Thus, 
we first set d = and solve for A^, and then use this 
value of As to solve for d. Setting d = in the above 
action, varying with respect to lA^p, and linearizing in 
small I As I by working near Tc, we obtain the linearized 
gap equation, 



kdk{ujl + k^ + H^) _ 1 

(^2 + fc2 + ^2)2 _ 4^2 fc2 - ^' 



where cj„ are fermonic Matsubara frequencies a;„ = 
(2n + l)TrTc. In the strong coupling limit g ^ 00 the 
presence of a non-zero fi is inessential and we can just 
use the strong coupling solution at /x = 0. However, in 
the weak coupling limit g — > the chemical potential is 
vitally important. At /i = 0, there is no superconductiv- 
ity, but as we shall discover, superconductivity develops 
for any non-zero /i. To see this, we sum over Matsubara 
frequencies to obtain the new gap equation. 



kdk 



tanh(|/c-|-^|/2rc) tanh (|fc - /x|/2Tc 



\k + fi\ 



7T 
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which defines Tc- Assuming fi <^ ujjj, and working in the 
weak coupling regime g/x <C 1, this has solution 



exp 



TT 



(13) 



Thus, we find there is a weak coupling instability to 
spin singlet superconductivity at non-zero temperature 
for any non-zero /i. However, the dependence on /i is 
strongly non-analytic. The expression ( 13 ) can be under- 



stood by noting that it is just the solution of the standard 
BCS Hamiltonian for a system with a low energy density 
of states proportional to ^. The zero temperature gap 
I As (T = 0)1 is proportional to T^. 

Thus, we have determined the magnitude of the pri- 
mary superconducting order parameter Aj. We now de- 
termine what happens to the triplet fields d in the pres- 
ence of the non-zero |As|. To see this, we expand the 



Lagrangian (12 1 in small d. Close to Tc (13), the expan- 



sion takes the form 



L = Lo+^a.y<dj +TrF4fcy^(A:d2(k) + Asd;(k)) 

- TrFAk,^i{Aldi{k) + Asd*i(k)), 

where Lq is the Lagrangian at d = 0, the matrix Uij has 
positive definite eigenvalues '-^ l/<?', and _F is a strictly 
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positive function of frequencies and momenta, which is 
even in w and k. The positivity of the eigenvalues of 
aij ensures that the critical temperature for spin triplet 
order to develop in the absence of spin singlet order is 
lower than the Tc (13). This condition must be satisfied 
in order for Ag to be the primary instability. 

Now note that at finite a non-zero A^ automatically 
generates a nonzero di^2 ^ fJ-^s- This can be straight- 
forwardly verified by minimizing with respect to di. In 
effect, As acts as a symmetry breaking field for the di at 
non-zero chemical potential. Moreover, the structure of 
the Landau expansion picks out di ~ k^/k = cos 0k and 
d2 ^ ky/k — sin^k, so that the integrals when multiplied 
by kx and ky respectively are non vanishing. Thus, we 
find that at small /x, the superconducting order parame- 
ter takes the form 



A+ =^(1'J(l + ii:5V(cos0k(7i -sincr2))iCT2^'_k), 

(14) 

where K is some undetermined positive constant and 
where g' is a function of f/,/i and ujd. The attraction 
in the triplet channel, g' , may already be non-zero at the 
bare level, and it can be verified that a g' interaction is 
also generated at second order in perturbation theory in 
small g. If = 0) = 0, then the leading contribu- 
tion to g' appears to scale as g' ^ [P'g^ juj^,. However, 
we have not verified this result and a rigorous determi- 
nation of g'{g^ fJ., ojd) lies beyond the scope of this work. 
Moreover, the behavior of g'{g, fi^ujjj) at strong coupling 
and near the critical point g = gc lies beyond the reach 
of perturbation theory in g, and may be a interesting 
topic for further work. Near the critical point g = gc, 
the expression (14 1, which is based on mean field theory. 



may also be invalid, and the nature of the singlet-triplet 
coupling near the critical point is an interesting topic for 
future study. For our present purposes, it is sufficient 
to note that we can smoothly go from the time rever- 
sal symmetric, gapped superconducting state at /i = 
to the time reversal symmetric, gapped state at large /i, 
according to ( 14 1. 



The phase diagram follows straightforwardly. The line 
segment ^ = 0, g < gc is semi metallic, and everywhere 
else in the n — g plane there is a single superconduct- 
ing phase, which is gapped and time reversal symmetric. 
The structure of the order parameter evolves smoothly 
according to 



^+ =H^*k(l+-F'(M)(cOS0k^Tl -sinCT2))iCT2^'-k>, 
k 

(15) 

where F(fi) interpolates smoothly between the limits 
F{^) = Kg' {g, ^,uj£i)ii for small ^, and -F(/i) = sign(/i) 
for large > ujd- 



sc 



SM 



gc 



SC 



FIG. 1: Phase diagram of clean Dirac fermion system. Along 
the line segment jj, — 0, g < the semimetal phase (SM) is 
stable. Everywhere else the ground state is a superconductor 
(SC) which has a fully gapped Bogolioubov de Gennes spec- 
trum and which preserves time reversal symmetry. There is a 
single superconducting phase, with an order parameter that 
is given by ( 15 1. 



III. 



MEAN FIELD SUPERCONDUCTIVITY OF 
DISORDERED DIRAC FERMIONS 



In this section we discuss the superconductivity of 
Dirac fermions in the presence of disorder. The section 
is structured as follows: first, we discuss disorder in the 
strength of the attraction g. We discuss the Harris crite- 
rion for determining the relevance of disorder, and show 
that the disorder is Harris relevant at the critical point 
/i = and g = gc- Thus, any disorder in the attraction 
strength must change the universality class of the crit- 
ical point. Next, we consider disorder in the chemical 
potential fi. We show that disorder in the chemical po- 
tential is relevant for g < gc on the fj. = line, and that 
it destroys the semimetal phase (Fig{T]) by introducing 
a non-zero density of states. The non-vanishing density 
of states in turn introduces a weak coupling instability 
to superconductivity. In the presence of chemical po- 
tential disorder, the zero temperature phase diagram in 
the entire n — g plane thus contains a single phase - a 
superconductor. We are thus driven to the remarkable 
conclusion that at fi — and g < gc, the presence of 
disorder actually enhances superconductivity. By intro- 
ducing a non-vanishing density of states, disorder allows 
superconductivity to develop in what would have been a 
semimetallic phase had disorder not intervened. 

The rest of this section is focused on understanding 
the behavior along the fj, = line, in the weak cou- 
pling regime, in the presence of disorder. We show that 
for smooth chemical potential disorder, where the disor- 
der is correlated over lengthscalcs large compared to the 
superconducting coherence length, the superconducting 
critical temperature Tc may be extracted from the clean 
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system results by treating the system as being 'locally 
doped.' Meanwhile, for short range correlated chemi- 
cal potential disorder (which we model using 'delta func- 
tion impurities'), we derive an estimate for Tc based on 
a 'uniform mean field' calculation similar in spirit to the 
standard Abrikosov-Gorkov theorj'^ for disordered super- 
conductors. This estimate for Tc is doubly exponentially 
small in the disorder strength. 



A. Disorder in g 

First, we consider static disorder in g i.e. we allow g to 
be spatially non-uniform, fluctuating about some mean 
value {g). The fluctuations in g are assumed to be in- 
dependent of time. The only critical point in the phase 
diagram FigJT] that is tuned by g is the critical point 
on the /i = line at g = gc, which was argued to dis- 
play emergent supersymmetry in the clean system.^l^ 
Whether disorder changes the universality class depends 
on whether the Harris criterion is satisfied.!^ The Har- 
ris criterion states that disorder changes the universality 
class if 



vd < 2, 



(16) 



where d is the spatial dimension and v is the critical 
exponent for the correlation length ^, which diverges near 
the critical point as ^ ^ (g — gc)~'^- Intuitively, if the 
Harris criterion is satisfied, then the typical fluctuation 
in g, averaged over a box of size is greater than the 
remaining distance to the critical point. 

For the particular critical point under consideration 
here, i/ = 3/4 (to leading order in an e-expansion) and 
d ~ The Harris criterion is therefore satisfied, and 
disorder in the interaction strength is a relevant pertur- 
bation. It thus follows that disorder in g should change 
the universality class of the interaction. A determina- 
tion of the 'true' critical point is beyond the scope of this 
paper. 



1. Smooth disorder 

The first model we will consider is of smooth disorder. 
Smooth disorder may be modeled by taking the Hamil- 
tonian Q at /i = and adding a term Hd, where 



Hd = J d^r^ir'iplao'ip: 

(Mr) = 0, 



r r 

2i?2 



(17) 



(18) 



The correlation length for the disorder R is assumed to be 
much longer than the superconducting coherence length 
specified in ( 20 ). In this limit, the system can be thought 
of as 'locally doped,' and consisting of a network of large 
electron and hole puddles. In each puddle, superconduc- 
tivity develops as if the system were doped with chemical 
potential ^o- In the weak coupling limit, the local order 
parameter at zero temperature is 



Ao LOD exp(-7rw^/5/io). 



which sets a coherence length 



V 

A 



■exp{Trv'^/gHo) 



(19) 



(20) 



This coherence length must be much smaller than the 
puddle size R to be in the smooth disorder regime. 

Thus, smooth potential disorder destroys the 
semimetal phase, replacing it by a network of large 
electron and hole doped puddles, each of which in- 
dividually becomes superconducting, with a critical 
temperature of order ([l9|. However, for the sample 
to be globally superconducting, it is essential that 
the various electron and hole doped puddles be phase 
coherent. 



2. Phase coherence of locally superconducting puddles 



B. Disorder in /i 

Static disorder in has an even more dramatic ef- 
fect: it destroys the semimetal phase by producing a 
non-vanishing density of states. The Harris criterion is 
the wrong framework for analyzing the effect of chem- 
ical potential disorder. The Harris criterion applies at 
a critical point which separates regions that are in dis- 
tinct phases. Meanwhile, the semimetal phase along the 
/i = line separates two regions which are in the same 
phase. As a result, all fluctuations about the semimetal 
phase place us in the same superconducting phase. The 
relevance of chemical potential disorder was also estab- 
lished using renormalization group arguments in Ref. 1121 
In this paper we will consider two distinct models of dis- 
order: smooth disorder and delta function disorder. We 
discuss each in turn. 



To estimate the critical temperature for phase coher- 
ence of locally superconducting puddles, we consider a 
specific model for disorder, which takes the form 



(21) 



/.tr = Ha sm{x/R) sin(y/i?). 



where R is much greater than the coherence length ([20j). 
This has electron and hole doped regions, as well as in- 
termediate regions which are close to undoped. The local 
coherence length scales as 



SO 



(ujd/v) 



where is defined by (20). Now the region near a; = 
with X < ^(x) can be thought of as an undoped 'barrier' 
region separating electron and hole doped islands. From 
this we conclude that the undoped barrier regions have 
width R/ln{R/^o). 
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Crucially, electron and hole doped puddles enter the 
same superconducting phase (Figjl]), which is mostly spin 
singlet A ^ {ipia2'4') for weak disorder /ig <C l^Jd- The 
Josephson coupling between regions of size R x R and 
local order parameter Aq, separated by a barrier region 
of width W iPl 



J = A, 



(22) 



Global phase coherence survives up to temperatures of 
order T fa J. Substituting W sa R/\n{R/£_o) into the 
above equation, we obtain an estimate for the critical 
temperature for global phase coherence, 



« Ao 



goln'^(i?/go) 
R 



(23) 



where Aq is defined by (19). Thus, global superconduc- 



tivity is weaker than local superconductivity by the small 
parameter £^o/R<^ 1. 

Note that the global Tc increases as R is made smaller, 
and appears to diverge as i? — > 0. However, the analysis 
is only valid in the smooth disorder regime ^o/R <C 1 
(and also in the weak disorder regime /io ^ i^d)- Thus, 
the global Tc is always smaller than the local puddle Tc 
by the small parameter £^q/R 1. 



3. Delta- function disorder 

In the strict weak coupling limit, g — 0, the coherence 
length (20) diverges, and disorder cannot be modeled as 



being smooth. There is therefore a need for a theory of 
superconductivity in the presence of short range disorder. 
The simplest possibility is to consider the limit i? -> in 

f. This may be modeled by adding to the Hamiltonian 
randomly placed positive delta function impurities 
with impurity strength V , and an equal number of ran- 
domly placed negative delta function impurities. 



H = Hi 



BCS 



(24) 



where rji = ±1, and V is the strength of the disorder. 

We assume that the 'impurities' are placed at random, 
so that the disorder concentration in a box of size L is 
Poisson distributed. 



1 



^/2TTno/L- 



: exp 



{n± - no/2)- 
no 



(25) 



Here PL{n±) is the probability that an i x L square 
box contains ± impurities with a concentration n±, uq 
is the mean total density of impurities. L must be much 
bigger than the typical distance between impurities, Iq — 

l/^A^o• 



C. Density of states from short range disorder 

In this section we discuss the density of states arising 
from short range disorder, in the absence of any interac- 
tion ((7 = 0). The density of states depends strongly on 
energy. For the clean system, the density of states scales 



^clcan 



27rw2 



(26) 



Meanwhile, an energy scale e also sets a wavelength 
A = vje. An electron with wavelength A 'probes' all 
impurities within a box of size A x A, and sees a local 
imbalance 8n = n+ — rt_ ^ 0, which give rise to a local 
chemical potential \i — VSn This in turn enhances the 
density of states according to 



V\Sn\ 



Now the probability distribution for 5n is 



P{5n) 



dnP{n-^- ~ n + (5n/2)P(n_ = n — 6n/2) 



exp 



2no 



which can be turned into a probability distribution for 
the correction to the density of states coming from a local 
imbalance 



XAtt 



V \f2jrriQ 



exp 



A 



2(2^^2)25^ 
1/22^0 



(27) 



The root mean square density of states from local imbal- 
ance scales as 



imbala 



A27rw2 ■ 



Thus, at the energy scale e — w/A, the correction to 
the density of states from local shifts of the chemical 
potential is 



imbalance 



27TV^ 



This scales with energy in the same way as (26). Thus, 



the 'local chemical potential' merely changes the co- 
efficient in (26) and does not qualitatively alter the en- 



ergy denendence of the density of states. However, there 
is a second contribution to the density of states, which 
comes about due to scattering of the electrons. 



1. Density of states from scattering 

The density of states can be extracted from the elec- 
tron Green function according to 



;y(w) 



1 



TrImG(a;), 
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where G is the retarded Green function. Now, the Green 
function for a Dirac fermion at ^ = 0, after ensemble 
averaging over disorder takes the form 



(k|G(c.)|k) = 



1 



where the scattering time r may be estimated using the 
self consistent Born approximation (SCB A) as in Ref. 1301 
This gives rise to a density of states equal to 



(28) 



where A is an ultraviolet cutoff of order the electronic 
bandwidth and no is the mean density of impurities. Note 
that there is a non-vanishing density of states even at 
zero energy. Thus, scattering on disorder destroys the 
semimetal phase, in agreement with Ref. [T^l The SCBA 
applies in the weak disorder limit noV'^ /v^ ^ 1. In the 
strong disorder limit the zero energy density of states 
scales linearly with impurity concentration, and can be 
thought of as coming from low energy bound states. In 
this paper we will focus on the weak disorder limit. 



D. Superconductivity from short range disorder 

We observed that in the presence of short range dis- 
order, the system develops a non-vanishing density of 
states. In the presence of a non-vanishing density of 
states, there is a weak coupling instability to super- 
conductivity. In this section, we search for a spatially 
uniform superconducting phase in the weakly disordered 
Dirac fermion system. 

We have a disordered system of Dirac fermions. The 
exact single particle eigenstates of the disordered system 
are \ipa)- The disorder is time reversal preserving. Thus, 
Kramer's theorem holds and all states come in Kramer's 
doublets \'4!a,a), where u is a Kramer's index. Let r be 
Pauli matrices acting in the space of the Kramer's dou- 
blet. The BCS Hamiltonian can then be written as 

Q 

where repeated indices are summed over. As usual, we 
have projected the interaction onto the BCS channel. We 
introduce a pairing field 

A/3 = ^gi3a{lpli<J2^a), 



by going to the path integral and using a Hubbard 
Stratanovich transformation. After decoupling the four 
fermion interaction, we obtain a Lagrangian of the form 



al3 



where a;„ is a fermionic Matsubara frequency. Upon go- 
ing to the Nambu spinor basis and integrating out the 
fermions exactly, we obtain an action purely for the or- 
der parameter fields, which takes the form 



L = -Trln[w2+e2 +|A„|2] 



E 

a/3 



From this we obtain the gap equation, which after inte- 
gration over Matsubara frequencies (at zero temperature) 
takes the form 



E 



A. 



/3 



(29) 



Note that we have not made any approximations in de- 
riving (29) (except for projecting the interaction on the 
BCS channel). So far, everything is exact, for a given 
realization of disorder. 

Now we ensemble average over disorder. After ensem- 
ble averaging over disorder, translation invariance is re- 
stored, and the eigenstates are indexed by momentum. 
The interaction is a constant in momentum space, so that 
the gap equation takes the form 



i'{e)de 



(30) 



where v is the disorder averaged density of states and g 
is the disorder averaged interaction. This 'disorder av- 
eraged gap equation' is the naive Dirac fermion analog 
of the 'Abrikosov-Gorkov' theory for superconductivity 
in disordered metals.!^ Now, the vertex correction to g 
arising from disorder at the one loop level is convergent, 
so that disorder does not produce a singular renormal- 
ization of g. We have g = Ag, where A is some 0(1) 
prefactor. We drop this prefactor for compactness, and 
use g — g- Substituting (281 into (30) and solving, we 
obtain 



« A 



uj D exp 



Lo £, exp 



gv 



exp 



(31) 



where the scattering time r is exponentially sensitive to 
disorder strength. Note that there is a non-zero for 
any value of however small. However, the critical 
temperature is doubly exponentially small in the disor- 
der strength. As we shall see later in the paper, (31 1 is 



a gross underestimate of the strength of the supercon- 
ducting instability in the disordered system. The true 
Tc is actually only exponentially small in weak disorder, 
not doubly exponentially small. However, we stress that 
Tc is non-zero even in (31) i.e. even though the clean 



system does not superconduct, the introduction of weak 
short range disorder introduces a weak coupling instabil- 
ity to superconductivity. This, disorder has the surpris- 
ing effect of enabling superconductivity, by destroying 
the semimetal phase. 
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IV. SUPERCONDUCTIVITY FROM RARE 
PUDDLES 

In this section, we investigate the possibihty that su- 
perconductivity may actually develop at a temperature 



much higher than (31 1, because of mesoscopic fluctuation 



eS'ects that are ignored in the uniform mean field analy- 
sis. The analysis in this section is inspired by the worlPH^ 
on superconductor to metal transitions in the presence of 
disorder, but with important differences arising from the 
different nature of the order parameter, and the fact that 
the 'critical point' /x = now separates two regions in the 
same phase rather than two regions in different phases. 

The specific possibility that we investigate is the fol- 
lowing: in a sample where disorder is weak and the den- 
sity of states is small, there may nonetheless be regions 
where disorder is stronger, and the local density of states 
is larger. These regions will have strongly enhanced local 
superconductivity. Even small fiuctuations in the disor- 
der concentration will have large effects on the local Tc, 
because of the double exponential sensitivity of Tc to dis- 
order concentration (31). Thus, one expects that in any 



disordered sample, there will be regions where the local 
Tc is much higher than (31). If the Josephson coupling 



between these locally superconducting regions is strong 
enough to establish global phase coherence, then the sam- 
ple will superconduct. The critical temperature for this 
'puddle based superconductivity' will be the temperature 
at which the locally superconducting regions lose phase 
coherence. As we will show, this temperature is only 
exponentially sensitive to disorder concentration, unlike 
the double exponential dependence in (|31[ ). Thus, it rep- 
resents a strong enhancement over ( |31[ ). 

This section is structured as follows: first we discuss 
the probability that a region of size L is locally supercon- 
ducting, and find the most likely size of superconducting 
puddle L*. Then we discuss the Josephson coupling be- 
tween superconducting regions, and the temperature at 
which global phase coherence is established. The discus- 
sion assumes that disorder is weak and dilute i.e. the 
sample is 'almost clean.' We conclude by making a few 
comments on the 'gauge glass' behavior that arises when 
external magnetic field is applied. 



A. Optimum size of superconducting puddles 

What is the probability that a puddle of a given size L 
displays local superconductivity at a given temperature 
T? Maximising this probability by varying L will tell us 
the optimal size of the superconducting puddle. 

We assume that the following equation is true 

Tc{y) w ud exp ( - \/gv{y)) , 

where T'c(r) is the local Tc in a region of size of order 
the coherence length, and v{t) is the local density of 
states in this region. A puddle of size L will be lo- 
cally superconducting if and only if the following two 




FIG. 2; Schematic illustration of puddle based superconduc- 
tivity. There are rare puddles which develop superconductiv- 
ity, and Josephson coupling then establishes phase coherence 
between these puddles 



conditions are satisfied: (i) The local Tc > T and (ii) 
the puddle is bigger than the local coherence length, 
L > (JIvf/i^d) exp(l/5i/). 

We can reformulate this as a condition on the density 
of states of the puddle. A puddle of size L will be locally 
superconducting IFF the density of states per unit area 
on the puddle 



> i^min — max 



1 



1 



ylnojjj/T' gln{ujDL/hvF 



(32) 



The probability that the puddle is superconducting is 
given by J^^,^ dvP{v). 

It is useful to define the thermal length 



Lt = v/T. 



(33) 



For puddles bigger than the thermal length, the temper- 
ature is the key cutoff, 

/•CO 

Psc{L>Lt) = / diyP{i^). (34) 

Jl/{g\nLT/a) 

Meanwhile, for puddles smaller than the thermal length, 
the puddle size is the key cutoff, 

/•OO 

Psc{L<Lt)^ / dj^P{iy). (35) 

Jl/{g\nL/a) 

Now, the density of states arises from disorder, and we 
can re-express 



PsciL) 



dvP{u) 



dnP{n), (36) 



where i/min = ^ max(j^^-^^, jj^^), and v is related to 
the local disorder concentration n by Eq. (28). In the 
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weak disorder regime, Eq. (28) tells us that lymin sets 

a.2 



V^lii{A/noV^iy„ 
Meanwhile, P{n) is given by 



1 



exp — L' 



no 



(37) 



(38) 



It is easier to integrate over P{n) since this distribution is 
just Gaussian. Upon substitution into (35 1 we find that 
the probability of superconductivity is 



Psc(L) = -erfc — 

2 V 



where we assume n,„in S> uq otherwise we are just dealing 
with uniform mean field superconductivity. Now nmin is 
given by ( 37 1 and Vmin is given by ( 32 ) . Thus, rimin has 
a log(log) dependence on L. Substituting into the above 
equation and plotting, we find that Pl{SC) is a monoton- 
ically decreasing function of puddle size over the entire 
range of sizes satisfying ngL^ 3> 1, with smaller puddles 
being exponentially more likely to be superconducting. 

Thus, smaller sized puddles are much more likely to 
superconduct. The smallest sized puddle that it is mean- 
ingful to talk about is a puddle with size of order Iq ~ 
1/y/no, where is the typical spacing between impuri- 
ties (which is assumed to be large in the weak disorder 
regime). Thus, the most probable superconducting pud- 
dle has a size of order Iq. The local Tc for this puddle 
may be determined by remembering that this puddle has 
a local coherence length (at zero temperature) of order 
Iq. Thus u/A sa v/Tc « Iq. This gives rise to a local 
critical temperature 



T, 



loc 



(39) 



In the weak disorder regime np — >■ this is small, but it is 
only linearly small in weak disorder, not doubly exponen- 
tially small. Thus, the local for the superconducting 
puddles is enormously enhanced over the uniform mean 
field Tc. 

The probability that a given region of size Iq supercon- 
ducts is 



PsM 



erfc 



erfc 



no) 



(no ln( 



noV^ 



- 1 



Now v'^/noV'^ S> 1 to be in the strong disorder regime, 
so the above probability is much less than one. We can 
approximate it as 



-Psc(^o) --exp - 



nlV^ In^ 



where we have neglected double log terms. The typical 
separation between superconducting puddles is then i?, 
where 



R ^ Iq exp 



1 



2 nlV^ In^ ( 



> ^0- 



(40) 



In order for the sample to be globally superconducting, 
the Josephson coupling between distant puddles must be 
strong enough to establish phase coherence. 



B. Josephson coupling between distant puddles 

In order for the puddles to be phase coherent, the ther- 
mal length must be larger than the typical inter-puddle 
spacing, otherwise thermal decoherence will destroy the 
Josephson coupling. This sets a bound on the tempera- 
ture for global phase coherence. 



KT 



< v^/no exp - 



1 



' exp 



(41) 



This is exponentially smaller than the local Tc ( 39 ) , but 



note that it is still only exponentially small in weak dis- 
order, not doubly exponentially small. 



As long as the bound (41 ) is satisfied, we can model the 



Josephson coupling between puddles using the zero tem- 
perature results obtained by Ref. [Ml In Ref. it was 
determined that the Josephson coupling between distant 
puddles at zero temperature scales as 



J 



where W Iq is the size of the puddle and R is the 
typical inter-puddle separation ( 40 ) , and the intervening 
region is modelled as being at the Dirac point. Note that 
the Josephson coupling is power law with distance, and 
hence long ranged. Taking into account the finite density 
of states in the intervening region will alter the power, 
making the Josephson coupling decay more slowly with 
distance, but we use the above result to be conservative. 
Taking into account the finite density of states in the in- 
tervening region will just add an 0(1) prefactor to the 
exponent in ( 42 ) (and will enhance the critical tempera- 
ture). 

Note that since the superconductivity is mostly s-wave 
close to the Dirac point, the Josephson coupling is un- 
frustrated. This is a major difference to Ref. HHSl where 
frustration arising from the d-wave nature of the order 
parameter dramatically impacted the physics. 

The unfrustrated Josephson coupling will be strong 
enough to establish phase coherence upto a temperature 
Tkt ~ J- This tells us that the maximum temperature 
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up to which global phase coherence can be expected to 



Vy/no exp 



In^ ( 



gv 



' exp 



gv ^ 



(42) 



This is a smaller temperature than the thermal bound 



(41 1, so we conclude that this is the true critical temper- 
ature at which global phase coherence is lost. Note that 
this temperature is only exponentially small in weak dis- 
order, not doubly exponentially small as in (31). Thus, 
the critical temperature arising from rare superconduct- 
ing puddles with phase coherence is enormously higher 
than the critical temperature for uniform superconduc- 
tivity. 



C. Magnetic field and gauge glass behavior 



(a) 



-■-o--- 



(b) 
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FIG. 3: One-loop diagrams in the theory with chemical po- 
tential disorder for (a) boson mass, field and velocity renor- 
malization and (b) fermion field and velocity renormahzation. 
The dashed lines represent boson propagators and the solid 
lines, fermion propagators. The couphng A is represented 
by four dashed Unes meeting at a point, and the disorder- 
induced coupling is represented by four soUd lines meeting 
at a point. A dashed line ending on a solid line corresponds 
to the boson- fermion coupling h. 



In the above discussion, the Josephson coupling be- 
tween superconducting puddles was unfrustrated, be- 
cause each puddle is mostly 's-wave.' The application 
of a transverse magnetic field frustrates the Josephson 
coupling, introducing a random phase difference / A • dr 
to each Josephson link, where A is the magnetic vector 
potential and the integral goes along the line connect- 
ing two puddles. At zero temperature, this turns the 
globa lly ph ase coherence superconductor into a 'gauge 
glass' .™lAt finite temperature, vortex creep (i.e., phase 
slips across the Josephson junctions) will introduce a non- 
zero resistance. Increasing temperature at non-zero mag- 
netic field then drives a smooth crossover to the high 
temperature semimetallic phase. 



(a) 
(b) 

-< = -ex 

FIG. 4: One-loop diagrams in the theory with chemical poten- 
tial disorder for renormahzation of (a) the 4-boson coupling 
A and (b) the boson-fermion coupling h. 



D. Fate of the triplet pairing mode 



In Section II, we showed that in the clean system away 
from charge neutrality /i 7^ 0, s-wave pairing induces spin 
triplet 'p-wave' pairing, with the p-wave amplitude being 
small when fj, <^ ujd- Insofar as a particular puddle has 
^ ^ 0, s-wave pairing will also induce some local p-wave 
component of the order parameter. However, this p-wave 
component will be small, and moreover the Josephson 
couplings between p-wave pieces will be strongly frus- 
trated, and so we do not expect any long range ordered 
p-wave component to the order parameter. 



perturbatively accessible in the e-expansiorP^I^ close to 
four spacetime dimensions. To take the disorder int o ac- 
count we use the replica trick.l^H We first show in Sec. V A 
using a one-loop RG analysis that chemical potential dis- 
order is perturbatively irrelevant at the supersymmetric 
critical point. However, this result is deceptive because 
chemical potential disorder generates disorder in the ECS 
coupling at the two-loop level, and this type of disorder 
turns out be relevant at the critical point as shown in 
Sec.lVB] 



V. REPLICA RENORMALIZATION GROUP 



In this section, we examine the interplay of disorder 
and superconductivity using a perturbative renormahza- 
tion group (RG) treatment. This technique is appropri- 
ate for studying what happens close to the supersymmet- 
ric critical point fi = 0, g = gc- The clean critical point is 



A. Chemical potential disorder 

We consider the problem of Dirac fermions ip coupled 
to the bosonic s-wave superconducting order parameter 



15 



with 



>■< * >--o--< • >--^----< 




FIG. 5: One-loop diagrams in the theory with chemical po- 
tential disorder for renormalization of the disorder strength 
A,. 



Ag. The Euclidean Lagrangian iJiUlSl 



4! 



where tp = ( ■0^ is a two-component spinor, c/ is 
the fermion velocity, Cb is the boson velocity, A > is 
a contact 4-boson interaction, /i is a boson-fermion cou- 
pling, and r ^ [gc — g) is a parameter which drives the 
quantum phase transition between the Dirac semimetal 
(r > 0) and the superconductor (r < 0). The Dirac 
conjugate is V' = — ^"0^70 where we choose the Dirac ma- 
trices to be 7o = fa, 7i = CTi, and 72 = 02, where cr^, 
i — 1,2,3 are the Pauli matrices. Short-range chemical 
potential disorder is modeled by a random variable /.<(x) 
with Gaussian distribution centered at zero. 



P[^(x)] cx e 



where ~ noV^ is the disorder strength. The disor- 
der can be integrated out using the replica trick."^ The 
replicated action takes the form 



n „ 
a=l 

a=l ^ 



n „ 



a=l ' 

A, 



X (0t^,)(x,r)(V'l^fc)(x,T'), (44) 

where n is the replica index to be set to zero at the end 
of the calculation. We denote the number of spacetime 
dimensions by Z? = 4 — e and the number of spatial di- 
mensions by d = _D — 1. The clean supersymmetric crit- 
ical point is accessible with a one-loop calculation,'2Sllll 
hence we perform a one-loop calculation including disor- 
der. The one-loop diagrams for the boson and fermion 
two-point functions are given in Fig. [31 and those for 
the four-boson vertex are given in Fig.TOa), the unique 
one-loop diagram for the boson-fermion vertex is given 
in Fig.[4][b), and the one-loop diagrams for the disorder- 
induced four-fermion vertex are given in Fig. [5] The lack 
of Lorentz invariance in the bare theory (c/ 7^ Cb) leads 
to anisotropic momentum integrals which are handled 
using an approach similar to that of Ref. 131] We write 
the spacetime momentum as (po, p) — where p is the 
magnitude p = ^/p^ + P'^ and n is a unit vector. The 
angular integral / (iO„ only contains information about 
the anisotropy Cf ^ Cb and does not diverge. We evalu- 
ate it in four dimensions, i.e., over the unit three-sphere. 
The integral over the magnitude p typically diverges in 
four dimensions and is evaluated in Z) = 4 — e dimen- 
sions. In Fig. [5] only the first three diagrams contribute. 
The remaining diagrams cannot generate an effective in- 
teraction which is nonlocal in time, and therefore do not 
renormalize A^. On the critical hypersurface r = 0, the 
one-loop RG equations are obtained as follows, 



dcj_ 
d£ 

dcb 
dX 



32h^{cb-Cf) 

3Cb{Cb + CfY 



S = Sf + Sb + Sbf + Sdis, 



(43) 



dh^ 



dA^ 

di 




5A2 192/1* 



Cb{Cb -f 

32/i2 



4/i^ 



Cb{cb + CfY 
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to 0{e). We search for the clean supersymmetric crit- 
ical point and linearize the RG equations around that 
fixed point. We indeed find a clean fixed point (A* = 0) 
with emergent Lorentz invariance (cj = — c = 1) and 
emergent N = 2 supersymmetry ((/i'^)* — j^, A* — e) 
corresponding to the Wess-Zumino model with one chi- 
ral multipletPS! A study of the full RG equations includ- 
ing the flow of r shows that the fixed point is indeed 
at r* — 0. Linearizing the RG equations at the critical 
point, we find one relevant direction, one marginal direc- 
tion, and four irrelevant directions. The relevant direc- 
tion corresponds to the mass parameter r with eigenvalue 
= 2 — e + 0{e^), which gives an order parameter ex- 
ponent v ~ i + I + ©(e^) in agreement with Ref. 231 
This is the direction which drives the transition. The 
difference of fermion and boson velocities Cf — Cb as well 
as the couplings A and /i^ have a nonzero projection only 
along irrelevant directions. (The sum of fermion and bo- 
son velocities Cf + Cb has a nonzero projection along the 
marginal direction.) The only direction along which the 
disorder strength has a nonzero projection is also ir- 
relevant, with eigenvalue 



= -1 



(45) 



which is negative and thus irrelevant for small e, and even 
in the limit e — >■ 1 corresponding to the physical case of 
2-1-1 dimensions (although corrections of 0{e^) cannot be 
neglected in this case). Therefore disorder in the chemi- 
cal potential appears to be an irrelevant perturbation at 
the supersymmetric critical point. 



(a) 



■o--- 



■-■0-- 



(b) 



FIG. 6: One-loop diagrams in the theory with boson mass 
disorder for (a) boson mass, field and velocity renormaliza- 
tion and (b) fermion field and velocity renormalization. The 
disorder-induced four-boson coupling Ay is represented by a 
square box. 
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B. Disorder in the BCS coupling 



FIG. 7: One-loop diagrams in the theory with boson mass 
disorder for renormalization of (a) the 4-boson coupling A 
and (b) the disorder strength Av- In this case there is no 
renormalization of the boson-fermion vertex at one loop. 



The analysis in the previous section could lead us to 
believe that the supersymmetric critical point is stable 
against disorder in the chemical potential. However, dis- 
order in the chemical potential will induce randomness in 
the coefficient of the |Asp term in the Landau- Ginzburg 
action ( 12 1 as can be seen from the presence of a /i^ | A^ p 
term in the clean case. In other words, an interaction of 
the form ( 44 ) but for the bosonic order parameter (f> will 
be generated at two loops, with a coefficient Ay c>c h^A^. 
This was missed in our one-loop calculation for chemical 
potential disorder. However, this interaction is pertur- 
batively relevant at the Gaussian fixed point and should 
be included in the calculation. More generally, random- 
ness in the BCS coupling g also gives rise to a random 
coefficient for the |Asp term. As a result, at the critical 
point, chemical potential disorder is a dangerous irrele- 
vant perturbation, which generates a relevant four-boson 
term. We therefore repeat the one-loop RG analysis but 
replace Eq. ( 44 1 by this four-boson term. The Euclidean 



action again takes the form (43), but with Sdis given by 



Ay 

2 



^ / ^'^'^ / / 



a,b=l 

X (0a0a)(x,T)(0^0fc)(x,T'), 

where Ay cx /i^A^ ~ n^V^ is the disorder strength. The 
one-loop diagrams for the boson and fermion two-point 
functions are given in Fig. [6j and those for the boson 
four-point functions are given in Fig. [7] Four-fermion 
interactions will be generated under the RG but are ir- 
relevant for small e. In contrast with the time-reversal 
symmetry breaking transition for Dirac fermions where 
the order parameter is in the particle-hole channel,'^^'^ 
here there is no renormalization of the boson-fermion ver- 
tex at one loop. On the critical hypersurface r = 0, the 
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one-loop RG equations are 

dcf _ 32ft,2(cb - cf) 
~dl ~ 

d£ 




to 0{e). We find the same clean (Ay = 0) supersym- 
metric fixed point as before. Linearizing the RG equa- 
tions around this fixed point, this time we find two rel- 
evant directions, one marginal direction, and three ir- 
relevant directions. One of the relevant directions corre- 
sponds to the mass parameter r with the same eigenvalue 
Ur = 2 — e + ©(e^) as before. The other relevant direc- 
tion is the only one with a nonzero component along the 
disorder strength Ay, and has the eigenvalue 



which is greater than zero and thus relevant for small e. 
In fact, all couplings (including the difference between 
boson and fermion velocities) have a nonzero projection 
onto this relevant direction. The criterion for the rel- 
evance of disorder i/Av > is equivalent to the Harris 
criterion, which in the context of the e-expansion should 
be written as d < 2yr with d = D — 1 = 3 — e and both 
sides of the inequality are expanded to 0{e). Although 
UAv appears to vanish if e is naively set to one to reach 
the physical case of two spatial dimensions, this is most 
likely only true at linear order in e, and for e — > 1 cor- 
rections of O(e^) and higher cannot be neglected. To the 
extent that the behavior at small e is representative of 
the physical problem, the signatures of the clean quan- 
tum critical point with emergent supersymmetry will be 
visible above a crossover temperature 



T* - AA 



1/(1- 
vo 



(46) 



where A is a high-energy cutoff which for the surface state 
of a topological insulator can be taken as the bulk en- 
ergy gap, and AyQ ^ tiqV^ is the bare disorder strength. 
For temperatures T < T* , the clean quantum critical 
behavior will be washed out by disorder and the dis- 
order Ay flows to strong coupling. Such a strong dis- 
order fixed point cannot be reliably studied within the 
present perturbative RG scheme. However, from a gen- 
eral standpoint we propose two possible scenarios. In 
the first scenario, the effective 4-fermion interaction in- 
duced by chemical potential disorder, which is irrelevant 
in 4 — e dimensions at the clean supersymmetric critical 



point [Eq. (|45|)], becomes relevant at the strong disor- 
der fixed point. In this case, we expect that a nonzero 
density of states would be generated for the fermions, 
and superconductivity would develop as a result of the 
Cooper instability below a nonzero critical temperature 
Tc < T* . In that sense, the crossover temperature (46) 



can be seen as an upper bound for at the critical point 
5 = (7c in the disordered system. 

In the second scenario, the chemical potential disorder 
remains irrelevant at the strong disorder fixed point. The 
density of states for the Dirac fermions remains zero at 
/J, = 0, and Tc is zero at g = gc- The behavior at g = gc 
would be controlled by the strong disorder fixed point. 
This would correspond to a non-monotonic dependence 
of the superconducting Tc on where Tc is zero at 5 = 
and g = gc, but nonzero for < 17 < ijc and for g > gc- 
Since such a non-monotonic dependence seems counter- 
intuitive, we expect this second scenario to be unlikely, 
and expect that the RG does eventually flow to a su- 
perconducting phase. However, we cannot exclude this 
possibility since the RG flows to strong coupling. 

We note that the replica field theoretic analysis as- 
sumes translation invariance and neglects mesoscopic 
fluctuation effects, which were shown to dominate the 
physics at weak coupling. If the true physics near the 
clean critical point is also dominated by mesoscopic fluc- 
tuations, then the replica field theory approach will dra- 
matically underestimate Tc- Conversely, if superconduc- 
tivity near the clean critical point is spatially uniform, 
and puddles are unimportant, then the nature of the 
superconductivity changes between weak coupling and 
strong coupling. This change in the nature of supercon- 
ductivity may then be controlled by a strong disorder 
fixed point. Similar ideas have been discussed for disor- 
dered bosons in one dimension in Ref l47l4Si We defer an 
investigation of these ideas to future work. 



VI. QUANTUM CRITICALITY IN THE 
DISORDERED SYSTEM 



It has been pointed out that the quantum criti- 
cal point in the ideally clean system is described by 
an unusual eff ective field theory displaying emergent 
supersymmetry-fiSlini However, we have also pointed out 
that the semimetal phase is itself unstable to disorder. 
Thus, disorder has the effect of destroying the quantum 
critical point discussed in Ref. I19I20I by inducing super- 
conductivity at arbitrarily weak interaction strengths. In 
this section we discuss to what extent it is possible to 
observe signatures of the (destroyed) quantum critical 
point. 

We assume that the system has been fine tuned to the 
Dirac point /i = 0. This assumption was also made in 
Ref. I19I20I The phase diagram of the clean system takes 
the form Fig.[8| Note the existence of a 'quantum critical 
regime' at finite temperature. In this regime, one can de- 
tect signatures of the proximate quantum critical point. 
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FIG. 8: Phase diagram of the clean Dirac fermion system at 
fi = 0. The system has only two phases - a superconductor 
and a semimetal, with a boundary between them which scales 
as Tc ~ (p — QcY" , where 2 = 1 from Lorentz invariance, and 
u = | + |+ 0(e^) in4 — e spacetime dimensions. Within the 
semi metallic phase there is a quantum critical regime, which 
is controlled by the critical point at g = Qc- The dashed 
line indicates a crossover between semi metallic and quantum 
critical behavior, and is given approximately by Eq. (|47|). 



The signatures of the quantum critical point gradually 
phase out as we move further from the quantum crit- 
ical point. The dashed line indicates the existence of a 
crossover between quantum critical and semi-metallic be- 
havior. Note that the quantum critical and semimetal- 
lic regions are not different phases. Rather the system 
evolves smoothly from one to the other. 

Where should one place the boundary of the quantum 
critical regime? We can answer this question for the clean 
system as follows. In general, the system may start any- 
where in the basin of attraction of the quantum critical 
point, but it will only start to display quantum criti- 
cal behavior when all the irrelevant couplings \i become 
smaller than some threshold small scale A*. The RG flow 
equations for the irrelevant couplings take the form 



~dl 



where the scaling dimensions yi were calculated in Sec. V, 
and / = InAo/A/. The RG is started at the initial scale 
Aq ~ LO]j (below which we have an attractive interaction), 
and stops at the scale Af Ri T, where the temperature 
T supplied the IR cutoff. From this, one obtains the 
criterion for quantum critical behavior. 



(47) 



where Ag is the bare coupling of the «th irrelevant opera- 
tor, yi is its scaling dimension, and A* is a small thresh- 
old. One should include only those irrelevant couplings 
which start far away from the critical point Aq > A*. 



FIG. 9: Phase diagram of the disordered Dirac system, with 
very weak disorder. The quantum critical point has been 
buried under the superconducting phase. However, signatures 
of the quantum critical point may still be visible in the finite 
temperature quantum critical regime. 



The precise choice of A* is somewhat arbitrary. For def- 
initeness, we suggest using A* — 0.1 as a threshold for 
quantum critical behavior. 

How does the phase diagram change in the presence of 
disorder? We have argued that in the presence of chem- 
ical potential disorder, the semimetal phase is unstable 
to superconductivity, with a critical temperature is given 
by either (23) or (42) depending on whether or not the 



disorder is smooth. Meanwhile, the critical point is also 
unstable to disorder, although the behavior at g = gc is 
controlled by a strong disorder fixed point which we were 
not able to access in any controlled manner. Although 
we cannot make definite predictions about the critical 
theory, since the RG fiows to strong disorder, we an- 
ticipate that the superconducting Tc should interpolate 
smoothly between weak and strong coupling. This leads 
to a phase diagram of the form shown in Fig|9] Note that 
the 'quantum critical point' has now been buried inside 
the superconducting phase. 

It is not possible to directly probe the quantum criti- 
cal point, since it has been buried by the superconduct- 
ing phase. However, at a temperature above the su- 
perconducting critical temperature Tkt, one can probe 
the quantum critical regime, to look for finite temper- 
ature signatures of the emergent supersymmetry. Since 
no real world sample is ever perfectly clean, it follows 
that signatures of the 'emergent supersymmetry' identi- 
fied in Ref. ll9l2Ul can only ever be probed by experiments 
conducted in this relatively high temperature quantum 
critical regime. 

We note that although the critical temperature for on- 
set of quantum critical physics is of order ujd (47), this 



boundary can be strongly suppressed if the bare theory 
starts a long way away from the critical point. In princi- 
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ical temperature is actually given by Eq. (42), which is 
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FIG. 10: Phase diagram of the disordered Dirac system, with 
less weak disorder. The quantum critical point and the quan- 
tum critical regime have both been buried under the super- 
conducting phase. 



pie, it is possible that the quantum critical regime may 
be entirely buried beneath the superconducting phase, 
in which case no signatures of the quantum critical point 
would be detectable in experiments. The resulting phase 
diagram would then look like (Fig. 



10). This scenario 



would arise if the critical temperature for quantum criti- 
cal behavior ( 47 ) were less than the critical temperature 



for superconductivity. 



VII. CONCLUSIONS 

We have calculated the phase diagram for a single 
species of Dirac fermions with attractive delta function 
interactions. We have shown that the zero temperature 
phase diagram in the coupling-doping {g — /i) plane con- 
sists of a single superconducting phase, except for a line 
segment along fi = and g < gc, which is a semimetal. 
We have shown that the introduction of disorder destroys 
the semimetal phase, introducing a finite density of states 
and triggering onset of superconductivity. Remarkably, 
disorder actually induces superconductivity, by introduc- 
ing a non vanishing density of states. This is a striking 
departure from the usual scenario, where disorder sup- 
presses superconductivity. We note that such a disorder- 
enhancement of superconductivity has already been ob- 
served in numerical simulations.^^ 

The critical temperature at which disorder enhanced 



superconductivity develops is given by Eq. ( 19 ) if the dis- 



order is smooth. The case of short range disorder is more 
complex. A naive estimate based on a uniform mean 
field solution for superconductivity gives an estimate 
(31), which is doubly exponentially small in weak dis- 



order. However, superconductivity is strongly enhanced 
by mesoscopic fluctuation effects, such that the true crit- 



only exponentially small in weak disorder. The super- 
conducting phase consists of locally superconducting is- 
lands, which establish global phase coherence through the 
Josephson coupling between them (Fig. [2| . Application 
of a sufficiently strong transverse magnetic field destroys 
the global phase coherence by frustrating the Josephson 
couplings, driving the system into a gauge glass phase. 

The region near the critical coupling g = gc was treated 
within a replica field theoretic approach. While this ap- 
proach neglects mesoscopic fiuctuations, it is useful for 
understanding the interplay of disorder and interactions 
in the strongly coupled theory. We find that chemical 
potential disorder is a dangerous irrelevant perturbation 
at the critical point, generating (relevant) disorder in the 
BCS coupling. This drives an RG flow to strong disor- 
der. While we cannot definitely answer what happens at 
strong disorder, continuity with the weak coupling results 
suggests that the RG flow ends up in a superconducting 
phase. As a result, the quantum critical point identi- 
fied in Ref. 1191201 will be buried under a superconducting 
phase for any non- vanishing value of disorder. However, 
for sufficiently weak disorder, signatures of the quantum 
critical point may still be visible in the flnite temperature 
quantum critical regime. 

This work also suggests some promising future direc- 
tions for research. Our analysis was focused on the prob- 
lem with a single Dirac fermion species, but graphene, 
a popular experimental material, actually possesses four 
species of Dirac fermions. A generalization of the anal- 
ysis to graphene would be a useful and worthwhile task. 
Another potential direction of research would be to fur- 
ther develop the analysis of mesoscopic fluctuations pre- 
sented in Sec. IV. While highly suggestive, this analysis 
was based on the assumption that the SCBA expression 
for the density of states could be applied at small length 
scales, to obtain a local density of states from a local con- 
centration of impurities, and that this density of states 
could be inserted into the BCS calculation. A more rig- 
orous treatment of this issue would be an important ad- 
dition to the present work. Meanwhile, the analysis at 
the critical point also opens up some avenues for further 
research. While we determined that the clean critical 
point is unstable to weak disorder, we found that the RG 
flowed to strong disorder. Determining what happens at 
strong disorder is a worthwhile topic for future work. In 
addition, the replica symmetric RG analysis ignored the 
effect of mesoscopic fluctuations, which were known to be 
important at the Gaussian point. Investigating the effect 
of mesoscopic fluctuations at the critical point is another 
topic for future work. Finally, an experimental investi- 
gation of the ideas outlined in this paper would present 
an excellent opportunity to compare theory with experi- 
ment. Given the rich phenomenology associated with su- 
perconductivity and disorder in Dirac fermion systems, 
and given the popularity of topological insulators as ex- 
perimental materials, we urge experimentalists to seach 
for superconducting topological insulators, and eagerly 
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anticipate further developments in this field. 

Note added: After completion of this work, we became 
aware of ReflSOl This work looked at the surface states of 
a topological superconductor with spin SU(2) symmetry, 
and concluded that they were unstable in the presence 
of vector potential disorder and interactions. It comple- 
ments our present work, which looks at the interplay of 
scalar potential disorder and interactions on the surface 
states of a topological insulator. 
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